Volume 99. Number 4, July-Augu&t 1994 

Journal of Research of the National Institute of Standards and Technology 

[J. Res, Natl. Inst. Stand. Technol. 9% 495 (1994)] 

Performance Comparison Between a 
Statistical Models a Deterministic Model, 
and an Artificial Neural Network 
Model for Predicting Damage From 
Pitting Corrosion 



Volume 99 



Number 4 



July-Augusi 1994 



M. Urquidi-Macdonald snd 
D. D. Macdonald 

Pennsylvania Slate Universiiy, 
University Park, PA 16802 



Various attempts have been made to 
develop models for predict ing the de- 
velopment of damage in mcinls and al- 
loys due lo pitting corrosion. These 
models may be di^^idc'd into two clov^es: 
the empirical approach which employs 
extreme value {Statistics, and the deter- 
ministic approach based on perceived 
mechanism!! for nuclcation and growth 
of damage. More recently, Artificial 
Neural Networks (ANNs), a nondeter- 
ministic t>T>e of motlel, has been devel- 
oped to describe the progression of 
damage due to pitting corrosion. Wc 
compare the three approaches above- 
statistical, dctcmiinistic, and ncunil net- 
works. Our goal h to illustrate the 



advantages and disadvantages of e^ich 
approach, in order that the most reli- 
able methods may be employed in fu- 
ture algorithms for predicting pitting 
damage functions for engineering siruc- 
turcs. To illusiraic ihc difficulty Ihat 
We face in predicting cumulative pitting 
damage, we selected a set of data ihat 
was collected in the laboratory. Wc 
«imparc and contrast the three ap- 
proaches by reference to this data set. 
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L Introduction 

On the basis of laboratory studies [1], and 
through the analysis of field data collected over the 
past decade by Batielle Columbus Laboratory [2], 
several factors have been identified as contributing 
to the development of pitting damage in gas fired 
heat exchangers in domestic and industrial service: 

(i) The type of alloy used for fabricating the 

heat exchanger 
(ii) Chloride concentration in the flue gas 

condensate 
(iii) Temperature 
(iv) Exposure time 
(v) Ambient versus indoor air 
(vi) pH 
(vii) Electrochemical potential 

Unfortunately, few of these factors are simply 
related to the damage functions or to one another. 



Accordingly, it is seldom possible to establish a 
.simple empirical equation for predicting pitting 
damage as a function of these variables. The case 
cited above is not atypical, and it illustrates the dif- 
ficulties laced by those who Seek to develop predic- 
tive models for assessing corrosion damage. 
Indeed, the data base established by Baltelle is 
probably one of the best that currently exists for 
the development of pitting damage in an industrial 
system. A full interpretation of the Batielle data in 
terms of statistical, deterministic, and artificial 
neural network models is published elsewhere [3]. 
In the present paper, we use a more restricted 
database to illustrate how various classes of models 
arc used to analyze the damage caused by pitting 
corrosion. These models include a statistical ap- 
proach based on the Weibull distribution function, 
a deterministic model based on a physicochemical 
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mechanism, and an Artificial Neural Network 
(ANN) that assumes neither a mathematical mode! 
nor a physical model, but which seeks to establish 
relationships between the dependent and indcpcn- 
dent variables by examining the patterns contained 
within the data set. 

2, Experimental Data 

We used laboratory data to illustrate the time 
and potential dependencies of pitting damage. To 
do this, we chose a laboratory data set for which the 
following independent variables were identified: 1) 
concentration of minor alloy elements in weight 
percentages, 2) difference in oxidation state be- 
tween host metal and minor alloy elements, 3) ap- 
plied potential, and 4) time of observation. Inde- 
pendent variables 1) and 2) are related to the type 
of alloy; and independent variable 4) together with 
solution composition (which was maintained con- 
stant) is determined by the electrochemical poten- 
tial. Temperature, solution composition, and pH 
were maintained constant The dependent variable 
was the total number of pits. 

We then used this set to illustrate the prediction 
of cumulative damage for pitting corrosion using 
three different models: statistical, deterministic, 
and artificial neural networks. The data were mea- 
sured by English and Macdonald at SRI Interna- 
tional [1]. 

Several alloys of nickel were fabricated. Each of 
the alloys tested was arc-melted from powders un- 
der an Argon gas blanket in a sealed container. Bi- 
nary nickel alloys containing Al, Ta, and Mo in 
nominal concentrations of 0.1 %, 0.5 %, 1 %, 3 %» 
5 %, and 8 % by weight were cast as 100 g buttons 
and were sectioned in an acrylic plastic before pol- 
ishing. The alloying elements were selected on the 
basis of their oxidation states relative to nickel (ox- 
idation state = 2). The excess oxidation states range 
from 1 for Al to 4 for Mo. 

The polished specimens were placed in a cell. 
The electrode potential was swept in the positive 
direction at 1 mV/s from an initial potential of 0.0 
V. This results in a distribution in breakdown po- 
tentials. Alternatively, the potential was stepped 
from 0.0 V to 0.325 V, 0.375 V, 0.4 V, and 0.45 V, 
This resulted in a distribution in induction (or ob- 
servation) times for the nucleation of pits. 

In both types of experiments the pit nucleation 
and growth events were photographed at 65 x mag- 
nification at regular intervals. The number of pits 
were counted on the pictures taken at different 
times -and conditions. 



The pitting data were measured several times on 
a similar sample to explore reproducibility. The re- 
producibility in pure nickel appeared satisfactory 
(about 10% difference between runs), but the 
reproducibility from alloy composition to alloy 
composition was different. Reproducibility was 
better at high potentials perhaps because the total 
number of pits developed was higher. Reproducibil- 
ity appeared to be better at high minor alloy con- 
tents and high oxidation states (about 20 %), than 
at low minor alloy contents and low oxidation states 
(about 50 %). Regardless of the poor reproducibil- 
ity in some of the samples, a general trend was ob- 
served: a) The cumulative number of pits 
diminishes with 1) minor alloy element content, and 
2) with increasing difference in oxidation stare be- 
tween the base alloy and the minor alloy element; 
and b) The cumulative number of pits grows with 
increasing applied potential and observation time. 

Cumulative phting damage is an irreversible, 
dynamic, time decay, environmentally related pro- 
cess. The literature is abundant in pitting corrosion 
data, but there is a lack of good quality data be- 
cause of the difficulty of measuring pitting corro- 
sion when controlling alt the environmental 
parameters. 

All model building is concerned with an attempt 
to increase our knowledge of complex physical real- 
ity. The parameters plus the validity of the model 
must be determined from the data. The philosophy 
behind the type of model is different. The informa- 
tion obtained from a purely probabilistic model 
(statistic and stochastic models) is about finding 
embodied in the data trends that can be used in fu- 
ture predictions. The information obtained from a 
deterministic model is about the physical meaning 
of the phenomena itself. The information obtained 
with a ANN model i.s about the dependency and im- 
portance of input /output relationships. In any case^ 
the model capabilities need to be tested. 

We can start the process of solving our problem 
by listing facts, listing observations, and listing exist- 
ing laws relating variables and outputs. Then we 
have to ask ourselves which will be the best model 
to describe the problem, and what do we expect 
frotn the model. Later we need to identify the 
model or models to use; specify the constraints, 
choose the coordinates, and apply the laws dictated 
by the model, important questions related to the 
choice of a correct model are: Is the process static 
or dynamic? Is the process stationary or not?; Are 
the available data distributed or not? What do the 
data mean? What is the data variance? What are 
the correlations?. In any case, the fitted model you 
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use for analyzing your data is the nearest icpicsca- 
tation of the true situation you have available. 

3. Statistical Approach and Results 

Stochastic processes are dynamic, and gix»d ex 
amples are fatigue, wear, and crack or pit groMili. 
There are two main types of stochasiic pro- 
cesses: stationary and nonstalionaiy. 

It is well known from experimental data thai cu- 
mulative pitting damage is a nonstationary phe- 
nomenon. It is well known that nojistationai^ 
models and their estimation are notoriously diffi- 
cult problems to handle except for special cases. 
Discrete stale continuous Markoff processes are 
good examples of models that describe nonsialiun- 
ary stochastic phenomena. However, there is no lit- 
erature on problem solving using nonsiationaiy 
finite Markoff chains [4]. On the other hand, for the 
lasl data set [1] (measured at the laboratory), ilie 
cumulative damage versus time was measured, but 
the pit depth versus number of pits was not. There- 
fore, it is impossible to derive a dynamic model for 
pitting damage using that data set. The only option 
available is to try to fit a static model (i.e., our hy 
pothesis is that the numbers of pits versus pit depth 
does not change with time). We choose a 2 p^iramc- 
ter WeibuH distribution; for which we assume thai 
ihe independent parameters are potential, and the 
oxidation state and concentration of the minor al- 
loying elements. The dependent variables are the 
cumulative number of pits and the induction time. 
The Weibull distribution function is 

F(A) = (l-exp(-(arW^)) 

where a and jS are fitting parameters and a- is the 
dependent variable of interest. 

We normalized the data set to 80 % of its maxi- 
mum value, allowing 20 % of the pits to nucleate if 
the time would have been extended to infinity. For 
each potential, oxidation slate, and percentile of 
minor alloy element, we performed a nonlinear fit 
to estimate the Weibull fining parameters. 

The choice of a Weibull distribution is arbitiary, 
we chose a Weibull distribution instead of some 
other probabilistic distribution because of Ihe flexi- 
bility that this distribution offers in fining dinercnl 
shapes obtained when plotting cumulative; damage 
versus dependent variables. 

The nonlinear fits were acceptable (Suu^ of 
square errors between fit and data <20 % for aov 
P). We used those data sets for which smooth 
changes of a and p were calculated as a fund ion ot 



puicniidl. iliLii 1 educcd the daia base to about 50% 
ut ihc loial axailahlc (the loial data base had 1400 
data lines coniaining number of piis al different ob- 
sorvaUvm linies, applied potential, oxidation states, 
and peiceiuilc uf minor elenjenis). We plotted the 
a-^uid jJ vahjes a^ d funciioji of potential. Figures la 
diid Ih 5.how the lesulis. The beta parameter of the 
Weibull dLsttihuiion appears to noi change with ap- 
plied pulcutial ai high concentraiion of minor alloy 
elements (5 %), but it changes drastically with ap- 
plied potential ai low concentrations of minor alloy 
elemcnl (3 %, 1 %), We then fit polynomials de- 
scribing a and fi as functions of applied potential, 
oxidatioij Slate, and percentile of minor alloy 
elements. 

"llie Wejbujl diiU ibuijon with a and J3 as parame- 
ters was used 10 generate the cumulative damage 
funcfion Figmes 2a, b, and c show the predictions 
obidined vvnh ihl^ siatisiical model. When we com- 
[iared the piediclions obtained with this model and 
the mezisurcd data, we observed that both trends 
are ^iniilar. lluwcvei", it would be very risky to use 
the model to mjke predictions for other oxidations 
states, percentile of minor alloys elements^ or 
applieJ purt-nlials outside the range for which the 
Weibull-paiametevs were calculated. 

It is well known that the Weibull distribution is a 
Siiffjoienily flexible function that practicaliy any set 
of dciia ccin be tilled by ji. However, the problem we 
fsLcd i^ thai uc do nut know u priori the correct re- 
lationship.-;; between a and fi and the independent 
variables. 

Picdicliuns wiih the same model for oxidation 
ataies greater than 3 2 gave cumulative probability 
of/Lio at any liinc and are not shown. The dcsigna- 
lii»n *'3 2" refcre.s to ihe oxidation state of the al- 
loying elenicnl (Al^.i) and the host metal (Ni=2). 



4, l)etcniiiui:^tic Model 

A Cuinplcicly sticcc.^.>rul model must account for 
all of the phenomenological correlations thai exist 
heiweeu pining susceptiVjility and pit velocity, and 
various tnviionnicnial and electrochemical factors, 
sueh as icmpeialure, pH, [CI], potential, time, and 
alloy Loniposiiion. The Poinl Defect Model (PDM) 
[5, 6] accounts for Ihe effects of electrochemical po- 
lejitial, alloy conii^ositioo, and environmental con- 
ditions 0!i llie nuclcalion of pits. 

The deterministic model is based on the PDM 
and the Solute Vacancy Inleraclion Model (SVIM) 
l7-i0]. The PDM propc^ses that passivity break- 
down occurs because of an enhanced flux of cation 
vacancies from the film/solution to the metal/film 
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Fig. 1, (a) Alpha parameler of the Wcibull distribution versus 
applied potential; at several percentiles of minor alloying cicmcntii, 
and oxidaliun stales, (b) Bctu parameter of ihc WcihuM di^itributton 
versus applied putentinl; at several perecntiles of minor alloying ele- 
ments, and uxidation States. 



interface. If the excess of vacancies arriving at the 
interface between the metal and the film can not be 
absorbed into the metal or be annihilated by some 
appropriate mechanism at high enough rate» they 
accumulate to form a vacancy condensate at the 
metal film/interface, which then grows to a critical 
size. The PDM is used to calculate the breakdown 
potential and induction lime. The effect of the mi- 
nor alloying elements in the oxide film on the 
breakdown parameters is modeled using the SVIM. 
The SVIM is based on the hypothesis that highly ox- 
idized solutes in the passive film electrostatically 
complex with the mobile cation vacancies. 

The PDM and SVIM results in distributed values 
of the breakdown potential and induction time, and 
complexing between the immobile alloying element 
in the film (the **solute") and the mobile vacancies 
diminishes the flux of vacancies across the film. This 
leads to an increase in the breakdown potential and 
the induction time for film breakdown. The higher 
the net oxidation state (minor alloy element oxida- 
tion-host ion oxidation) and/or the higher the per- 
centile of minor alloying elements in the film, the 



greater the effect on reducing the flux of vacancies 
and hence in increasing the pitting potential and 
the induction time. Once the pits nucleate, they 
grow at different rates. To calculate the pit growth 
rate we used 1) a simplistic steady state model sug- 
gested by Alkire [1 1]; and 2) a nonstationary model 
developed by us [12-13]. The stationary model is ex- 
pected to be adequate for only short times. 

The overall model (combination of the PDM, 
SVIM, and pit growth) requires the defining of a 
number of parameters, as shown in Table i. 

Figure 3a shows the cumulative probability of the 
number of pits (normalized to 1) as a function of pit 
depth and observational time of 50 s, for an applied 
voltage of 0325 V and for several concentrations of 
the minor alloying element with oxidation state of 
3-2 (example aluminum in nickel). Figures 3b and 
3c show similar plots for oxidation states of 4-2 
(e.g., titanium-nickel) and 6-2 (e.g., molybdenum- 
nickel), respectively, for the same conditions. It is 
interesting to note the great effect of minor alloying 
elements with high oxidation states. The model pre- 
dicts that the cumulative probability of the number 
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Fig, 2* Predicted cumubLivt: probabih'ty, ebtuincd using the statisti- 
cal model, versus lime of observation, at several applied poteniinls. (a) 
Oxidation state 3-2; and 1 % of minor alloying cicmcju segregated in 
the film, (b) Oxidation state 3-2; and 3 % of minor alloying element 
segregated in the film, (c) Oxidation state 3-2; and 5 % of minor 
alloying clement segregated in the film. 



Table 1. Input data used in the calculation of the deterministic/probabilistic model 



Parameters 



Value 



Units 



Stoichiomctry 
Avogrado constant 
Mol vol of oxide cation 
Gibbs energy change" 
Gibbs energy change* 
Mean diffusion coefficient 
Standard deviation 
Chloride activity 
Electrical field across film 
Alpha 
Beta 

Critical area vacancy size' 
Critical vacancy flux" 
Temperature 
Applied potential 
Molar gas con^^tanl 
Electrical pott^ntial film/sol 



6.023 E + 23 

3a 

-10»000 

5E-20 
0.75 Dmcan 
0,573/2 
1.1 E+6 

0.65 
-0.01 

1 E+16 

15.87 E+ 12 

298,15 

"0^5 

8.314 

-^0.5 



mol-* 
cmVmoi 

J/moI 
J/mol 
cm-/s 
crnVmoI 

V/cm 

V/pH unit 

No7cm- 

N'o7cm- ' ^ 

K 

VSHE 

JK-^mo|-* 

V SHE 



.■Variables that were used to adjust one datum point to scalt; the rcsuUi* properly. 
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Fig- 3a* Cumulative probability, calculated using the dclerministic 
model, versus pll dept^ for sc\'cral conccnlratiojis of minor allaying 
clement segregated in the film for oxidation state (3-2). 
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Fig. 3b* Cumulative prabability. culculatcd using the deterministic 
model, versus pit depth for several concentrations of minor alloying 
clement segregated in the film for oxidation state (4-2). 
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Fig. 3c. Cumulative probability^ calculitlttd using the deterministic 
models versus pit depth for several conccntraiions of minor alloying 
element segregated in the film for oxidation state (6-2). 



of pits at all pit depths is higher at lower minor 
alloying element oxidation state and at lower con- 
centration of the minor alloying element. Because 
the model does not assume a total number of break- 
down sites only a normalized probability is ob- 
tained. 

Figure 4 shows the beneficial effect of adding 
minor alloying elements with high oxidation states. 
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Fig. 4. Difference between cylculated breakdown potential 
of nickel containing %-5 % of minor alloying elements with 

oxidation slates of ( 3*2: Ni-A1), (---4-2: Ti-NI). and 

(-.-.-.6-2: Mo-Ni) and calculated breakdown potenitai of 
pure nickel (containing % of minor alloying cicmciits). 



The breakdown potential is shifted in the positive 
direction, indicating that higher potentials are nec- 
essary to achieve the same damage. 

S. Artificial Neural Network Model and 
Results 

Probably the most efficient method, when data 
are available, of establi.shing relationships between 
inputs and results is to use artificial intelligence 
techniques. Accordingly, we describe here an Artifi- 
cial Neural Network (ANN) for predicting pitting 
damage functions for condensing heat exchangers. 
When the net is trained with reliable data and 
knowledge, we are able to accurately predict dam- 
age outside the ranges of the input variables. 

An ANN is a highly interconnected system 
inspired by the brain and formed by simulated 
"neurons" represented by a transfer function, and 
"weights" associated to the connections of the 
"neurons." The back propagation training al- 
gorithm allows experimental acquisition of input/ 
output mapping knowledge within multilayer 
networks. Because wc have experimental data on 
the cumulative numbcrsof pits versus time of obser- 
vation, as a function of oxidation state, minor alloy- 
ing element, and applied potential, we decided to 
use an ANN backward propagation technique with 
supervised learning. During training of the ANN, 
the cumulative numbers of pits were used as 
'^output" and the applied potentia), oxidation state. 
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minor alloying element concentration, and time of 
observation as "inputs*" We explored several to- 
pologies to obtain the best compromise between 
learning and computing time for an ANN with 2 
hidden layers. 

The maximum training time was set to 12 hours 
on a Macintosh 11 microcomputer with a threshold 
of 10 % of the normalized input values (input- 
ouput) [2]. 

The ANN had the following features: 

(i) Heteroassociative memory, for which the 
patterns on recall from the memory are purposely 
different from the input pattern, because the inputs 
and outputs are different and belong to different 
classes of information. 

(ii) Delta rule type of learning, where the neu- 
ron weights are modified to reduce the difference 
between the desired output and the actual output of 
the processed element. The weights are changed in 
proportion to the error calculated. This rule also 
limits the learning, if the error at the output of the 
network is lower than a given threshold. The learn- 
ing rates of those layers close to the output are set 
lower than the learning rates of the other layers. 

(iii) A momentum term, which is used to 
smooth out the changes. 

(iv) A sigmoid transfer function, which is a 
monotonically continuous mapping function. 

The ANN predictions are in good agreement with 
the measured data. Figure 5 shows that correlation. 
Considering the difficulty of obtaining high quality 
data, we consider that the correlation factor is satis- 
factory. 



*5 




-4-2 5 4 6 

LN (Msasured Cumulative Mumber of Pas] 

Fig. 5. NdturdI lug^rilhm of the predicted ANN total number 

of pits versus natural logarithm of laboratory measured total 
number of pits. The measurements included several: applied 
potentials, observational times, oxidation states, and percent of 
minor alloying elements. 



After the ANN was trained, it was used to make 
predictions of the number of pits at different ap- 
plied potentials, observation times, oxidation states, 
and percentages of minor alloying elements in the 
film. The total number of pits predicted by the 
ANN decreased with increasing percentage of mi- 
nor alloying elements in the film, and with increas- 
ing oxidation state of the alloying element (Figs. 6a 
and 6b). Behavior similar to that predicted by ANN 
was observed experimentally. 

The ANN, once trained, can be used to explore 
the importance of the relationship between 
"output'' and "inputs." We found that the results 
were strongly dependent on observation time (/\ t), 
have a medium dependency on applied potential 
(K^'^), and show weak dependencies on oxidation 
state (Z'^) and concentration of minor alloying 
element in the film ([%]^''). 
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Fig. 6. ANN prediction of cumulative numbcrof pits versus ap- 
plied potential^ at 50 s ifme of observaljon; and several pcrcen- 
tiles of minor alloying elements. (6a) Oxidation state (4-2). (6b) 
Oxidation state (6-2). 
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6. Discussion 

Figures 7a to 7d show the best predictions 
obtained with the three models compared with the 
laboratory data. 

The deterministic model predicts that the cumu- 
lative probability at low applied potentials is not 
only described by a flat curve but that the curve is 
displaced to higher times. This prediction coincides 
with the experimental observations. Tht predictions 
with the deterministic model at high potentials indi- 
cated that the plateau corresponding to higher 
times is reached sooner than that measured. The 
deterministic model is the only model (compared to 
the other two models) that brings together an un- 
derstanding of the problem as well a predictive tooL 
Another advantage that the deterministic model 
has over any nondeterministic model is that to fit 
the model, only an experimental datum point is nec- 
essary to calibrate the model to the data. The deter- 
ministic model was developed to predict damage 



and cumulative number of pits simultaneously. This 
last capabih'ty makes it very attractive to the user* 

The results obtained with the probabilistic model 
are in general agreement with the experimental ob- 
servations. As with the deterministic model» the 
plateau in cumulative damage is reached sooner 
than the measured one. However* the curves are 
flatter at lower potentials than at high potentials, 
but they are not displaced to higher times. The 
probabilistic approach needs a large data base, and 
the predictive capabilities are limited to the ranges 
of variables confined in the data base. 

The ANN model describes the cumulative num- 
ber of pits very close to the experimental measure- 
ments. The plateaus on cumulative damage 
correspond very well to the plateau obtained exper- 
imentally. The ANN predictions at low number of 
pits are inaccurate, but they are very close to the 
experimental observation at higher cumulative 
numbers of pits. 
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Fig. 7* Cumulative number of pits versus time of observation at several jipplicd potentials; Oxidation ^tatc (3-2); and 3% of minor 
'dlloying clement. (7a) Predictions using the deterministic model. (7b) Predicdons using Wcibull Distribution model wi(h 2 filling 
paranficters, {7c) Prcdielions using (he ANN, (7d) Laboratory mcasuremenls. 
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We conchide that Utile is learned abnm xUf* phe- 
nomena when nflndetcrmini<;uV mod^H pif u'^erl; 
however, they can represent mvnlnMMe nr^t'liftive 
tools. The stalistfcal model is in a s<.m>';<' m^re de 
manding that an ANN rnodcK h re'ji'"'^*^ a uil»u^t 
and large data base. Wf^ found \hM mi A,|JN crn 
learn from "noisy" data, and thnt i^e ranee r*! j>ir 
diction can be extended outside the *nfige f^u which 
it was trained, if trained correctly \\}X Damnjge 
functions can also be calculated using i he deierniin- 
istic model based on ihe PDM and SVIM. Tlir^ *\e 
terministic model does not tieed to have ^n 
extensive data base that inHiides pit depth di*:tTihu 
tions. On the other hand, the nondrTrriulni<;iir 
models need a large data base, and (hey nre inKihle 
to make predictions of partial damage If »hc pit 
depth versus number of pits is not confined to the 
data base. Cumulative damage can be interpolfited 
and extrapolated to other voltages and rimes, and lo 
other applied potentials, for any of the thicc m'Kl- 
els* In general the results obtained with the. three 
models were found to be in reasf»ii:»ble njjfV'pniP'nt 
with experimental data [12]. 

We do not intend to emphasize here the impor 
tance of deterministic models over riondetepninls- 
tic models, but we have to keep in niiiui. when 
picking a model, to choose the. one thnt Hesi iTpre- 
sents the observations and that is reasounbly ra*;y lo 
solve. Clearly, the reliability of the p^inp^hitiori. ii> 
particular, depends critically on the f]n^liiv of tli-: 
data and on the veracity of the. model 



[7] D. D. Macdtinatd and M. Urquidi-Macdonald, Elcc- 

in^hrm. Acta 3t, UW (1986). 
r?^! M, llrt|riuli-Mricd*»ti:»ld and D. D. Macdonald» J. Elcc- 

irorhfin. Six:, 132 (3). S55 (1986). 

(o] M. tJrqiiidi-Macdonald and D. D. Macdonald. J. Etcc- 
trL^chcm. .Soc. 132, 5*55 (1^85). 
r n»l T,\ D, Mnotonaid, M. Urqnidi-Mycdonatd. and S. Lenhart, 

I P n L R F*cck a-id R. C Alkire, J. Electrochemicat Soc. 126, 

K>rj2(1979). 
['7] D. O. Macdunald and M. Urquidi-Macdonald, Corrosk>n 

46. 381) (1990). 
[^■^1 Nt. Urtiiiidi Macdonald, unpubllsticd data, Engineering 

S<*fcnco and Mrchimics^ Penniytvania State University 



About ihe authors: A/. Urquidi-Macdonald is an as- 
soriated professor at the department of Engineering 
Science and Mechanics of Fenn State. Her primary 
interpsis are in mathematical modeling of damage and 
life prediction. IX D, Macdonald is the Director of 
Fcnn State*s Center for Advanced Materials and is a 
Professor of Materials Science and Engineering. One 
of his principal research interests is the development of 
detninmisi'w models for predicting damage due to 
loni!i7ed connsinn. 



Acknowledgements 

The authors gratefully acknowledge the Mipp'Mt 
of this work by the Deparrmeni ryt i'ne*g\ 'Ft;i«ri'.' 
Energy Sciences. 



7. References 

[1] SRt Finiil Report. June I9fi9-Jiiry 1991. MntcriaT Rr^*>rif.-h 
Laboratory, SRI Intcmarlonnl, M^nl^ t'iirk. TA. 

[2] R. Razgaitis, J. H. Payer. S. G. Tnttirrt. fi. llmdiM. H \ 
White, D. W. IjiKklin, R. A. CudniV, and Q \\. .Ktikftml, 
Condensing Heat Exchanger for Rcsid^ntial/rotnmerci^d 
Furnaces Boiters, ?\\v\^ IV, P^itotlc R'-fKni tn vnyf 
BNL, BNL Report NVk 5t94 

[3] M. UrquidJ*Macdon.iPd And D. O. Mjciltmnld. ';**T(ii^inT\ 
to be submitted for ptibltcation f1*^94). 

[4] J. L. Bogdantiff and F. Kozin, Prcvb.'itiiiictic M^M--k fi 
Ciimnlative Damage, John Wilev A Soim ( IWO 

[5] C. Y. Chao, U F. [Jn, and n D M-.r.|.^,i^T.t. t Fji^r 
irochem. St>c, 128, 1187(1981} 

[6] L. F. Lin, C Y. Chad, wnd D P Nt^cd^v,-.td, T TK^. 
irochem. Soc. 12«, liui (t<},SO. 



504 



